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We study the quantum dynamics of a Bose Josephson junction(BJJ) made up of two coupled 
Bose-Einstein condensates. Apart from the usual ac Josephson oscillations, two different dynamical 
states of BJJ can be observed by tuning the inter-particle interaction strength, which are known as 
'•7r-oscillation' with relative phase it between the condensates and 'macroscopic self-trapped' (MST) 
state with finite number imbalance. By choosing appropiate intial state we study above dynam- 
ical branches quantum mechanically and compare with classical dynamics. The stability region 
of the '7r-oscillation' is separated from that of 'MST' state at a critical coupling strength. Also 
a significant change in the energy spectrum takes place above the critical coupling strength, and 
pairs of (quasi)-degenerate excited states appear. The original model of BJJ can be mapped on 
to a simple Hamiltonian describing quantum particle in an 'effective potential' with an effective 
Planck constant. Different dynamical states and degenerate excited states in the energy spectrum 
can be understood in this 'effective potential' approach. Also possible novel quantum phenomena 
like 'macroscopic quantum tunneling' (MQT) become evident from the simple picture of 'effective 
potential'. We study decay of metastable '7r-oscillation' by MQT through potential barrier. The 
doubly degenerate excited states in the energy spectrum are associated with the classically degen- 
erate MST states with equal and opposite number imbalance. We calculate the energy splitting 
between these quasi-degenerate excited states due to MQT of the condensate between classically 
degenerate MST states. 

PACS numbers: 03.75.-b,74.50.+r,05.30. Jp 
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I. INTRODUCTION 



The creation of Bose-Einstein condensate (BEC) of 
alkali atoms [l|, Q has opened up various possibilities 
to study coherence properties of macroscopic 'matter 
wave'. One of such direction is to study the Joseph- 
son oscillation^ in two weakly coupled BECs. Simi- 
lar to the superconducting Josephson junction, a Bosc- 
Josephson junction (BJJ) can be created by two BECs 
in a double well potential, which are weakly coupled 
by the overlap of their wave functions 0, [f|. The con- 
densates in each well are described by a 'matter wave' 
with a well defined phase. After they are weakly con- 
nected, the relative phase changes and a particle current 
is produced across the barrier. The nonlinear dynam- 
ics of ac and dc Josephson effects of BJJ have already 
been studied experimentally by loading condensates in 
a double well trap@, 0|- Also the matter wave inter- 
ferometry have been studied in a double well geometry 
on an atom chipQ. The number difference between two 
Josephson coupled condensates and their relative phase 
are two canonically conjugate collective variables which 
show coherent oscillations when an initial number imbal- 
ance is created between the condensates. The non-linear 
dynamics of the relative phase and number imbalance be- 
tween the coupled condensates have already been studied 
theoretically by several authors within the classical field 
approximation [9l-|l2| ■ Effect of dissipation in BJJ has 
been studied theoretically [iH lh"H and also been observed 
experimentally in the decay of Josephson oscillation [tJ. 



Quantum fluctuations become important for small num- 
ber of particles in the condensate and for large interparti- 
cle interaction. Full quantum dynamics of BJJ with finite 
number of particle shows various interesting non-linear 
effects like 'revival' and 'beating' for different values of 
coupling constant [TH, HU. 

One of the main advantages of BJJ is the possibility 
to tune the interparticle interaction strength which en- 
ables to study various types of non-linear Josephson os- 
cillations by controlling the interaction strength. Apart 
from the usual ac Josephson oscillation two other types 
of dynamical branches can exist and their stability re- 
gion is separated at a critical coupling strength. Above 
the critical coupling strength the 'macroscopic self- 
trapping' (MST) can occur dynamically with non vanish- 
ing time average of the population imbalance between 
the condensates [I]- Below the critical coupling strength 
another dynamical state known as V -oscillation' can ex- 
ist in which the relative phase between two condensates 

iS 7T |,[H|. 

We also notice a remarkable change in the energy spec- 
trum associated with the change in dynamical branches, 
and doubly degenerate excited states appear above the 
critical coupling. The system of many-bosons in BJJ can 
be mapped on to a quantum mechanical problem describ- 
ing a particle in an effective potential. The main focus 
of this work is to analyze different types of dynamical 
states and change in the energy spectrum in terms of 
this 'effective potential'. We also address novel macro- 
scopic quantum effects related to 'macroscopic quantum 
tunneling' (MQT) in this system which is evident from 
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the method of 'effective potential'. The decay rate of the 
metastable state corresponding to V -oscillation', and 
tunneling time between two MST states with equal and 
opposite number imbalance, are calculated analytically 
within semi-classical formalism of MQT. The MQT is 
always suppressed by the number of particles in the con- 
densate, nevertheless this novel quantum effects can be 
observed in B J J with small number of particles and close 
to the critical coupling strength. 

This paper is organized as follows. In section II, we 
model B JJ within single mode approximation by two site 
Bose-Hubbard model to study Josephson dynamics. By 
mapping the original Hamiltonian to a spin model we 
recover the 'pendulum model' @ describing the classi- 
cal dynamics of the relative phase and number imbal- 
ance. Different dynamical branches and collective fre- 
quencies of small oscillation, obtained from the classi- 
cal analysis are compared with the full quantum dynam- 
ics. The equivalence between the spin model of BJJ and 
the quantum mechanical problem describing a particle 
in an "effective potential" has been outlined in section 
III. Dynamical states below (and above) the critical cou- 
pling constant are analyzed from the effective potential. 
Next we outline possible novel collective quantum effects 
which emerge from the simple picture of 'effective po- 
tential'. In subsection A, we consider possible decay of 
meta-stable V -oscillation' due to MQT through the po- 
tential barrier. Subsection B is devoted to analyze the 
energy spectrum of the Hamiltonian and its connection to 
the self-trapping phenomenon. The appearance of dou- 
bly degenerate excited states in energy spectrum above 
the critical value of coupling can be understood in terms 
of a double well 'effective potential' whose minima corre- 
sponds to the dynamical fixed point of self-trapping with 
equal and opposite number imbalance. In subsection C 
we consider the MQT of the system between these two 
MST states and calculate the energy splitting between 
the quasi-degenerate excited states. Finally we summa- 
rize in section IV. 



II. THE MODEL AND JOSEPHSON DYNAMICS 

Bose Josephson junction of two weakly coupled BECs 
can be well described by a two-site Bose-Hubbard 
model(BHM), 

H = -J(a\a 2 +a\ai) + ^ [hi{h\ - 1) + h 2 {fi2 - 1)] (1) 

where aj(a,) is creation(annihilation) operator of bosons 
at site index i (i = 1, 2 for two sites), hi = aja, is number 
operator, J is the tunneling matrix element, and U is the 
on-site interaction strength. The tunneling matrix J and 
the on-site interaction strength U are the parameters of 
the effective model, which can be calculated from the 
exact shape of the trapping potential and macroscopic 
wave function of the condensates hfj. For a fixed value 



of total number of particles N in the condensates, we 
can transform the above Hamiltonian in Eq.([T]) to a spin 
model, 

H = -JS X + |s 2 , (2) 

where, J = 2J, U/2 = U, S x = \{a\a 2 + 4 a i); 

S y = ^(a\a2 — a\a\), and S x = ^{a\a\ — a^a^) are the 
Schwinger boson representation of a large spin of mag- 
nitude S = N/2. In the effective spin Hamiltonian we 
neglect the constant term depending on the total num- 
ber of particles N. 

For a large spin of magnitude S, we can 
represent it classically by a spin vector S = 
S(sin 9 cos 4>, sin 9 sin</>, cos 0), where two angles 9 
and <j> describes the orientation of the spin in the carte- 
sian coordinate system. Physically, the variables cos0 
and <j) describe the fractional number imbalance and 
relative phase between two condensates in the original 
BJJ model. In terms of these variables the Lagrangian 
of the classical spin is given by, 

L= S(cos9)4>- H(6,4>) (3) 

where <f> is the time derivative of </>, and H(9, (f>) — 
— JS sin 9 cos </> + cos 2 9 is the classical energy corre- 
sponding to the Hamiltonian in Eq.([T]). Here the variable 
cos 9 is canonically conjugate to the variable <f> and plays 
the role of momentum. To use dimensionlcss variable, we 
scale time by 1/J, energy by J and introduce a dimen- 
sionlcss coupling constant a = ^j-, in order to make the 
total energy an extensive thermodynamical quantity. 

From the Lagrangian given in Eq.Q, we obtain the 
classical equations of motion of a large spin, 

9 = sm<t>, (4) 
sin 6*0 = cos 9 cos cf> + a cos sin 0. (5) 

The linear stability analysis of the above dynamical equa- 
tions reveals stable dynamics around three fixed points, 
given below: 

• First fixed point is at 9 = it/ 2 and = 0. It de- 
scribes the ground state of the total system, which 
is the symmetric combination of the wave func- 
tions of two equally populated condensates in two 
well. The frequency u> of small amplitude oscilla- 
tions around this fixed point is, 

uj/J = \/TTa (6) 

which is the plasmon frequency of ac Josephson ef- 
fect. 

• A 'ir-phase oscillation! is described by the second 
fixed point at 9 = tt/2, and <f> = tt 0. This fixed 
point represents the wave function of first excited 
state which is the anti-symmetric combination of 
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wave functions of two condensates with equal num- 
ber of atoms and relative phase tt. The frequency 
of oscillation around this hxed point is given by 
lj/J = \/l — a. This '7r-oscillation' becomes dy- 
namically unstable when a > 1. 

A third fixed point at sm8 = — , and <f) = it 
represents a 'macroscopic sclf-trapped(MST)' state 
with non-zero population imbalance between two 
condensates Classically this self-trapping phe- 
nomenon can be understood from a Ginzburg- 
Landau(GL) type potential. A second order differ- 
ential equation for population imbalance q = cos 9 
can be obtained from the Eq.Qand Eq.©, 



dV GL (q) 



(7) 



V GL (q) = l(l-Ea/S)q 2 + ^q i (8) 

where q denotes second derivative of q with respect 
to time, the conserved classical energy is given by 
E/S = — yl — q 2 cos <fi + ^q 2 . Above equation for 
number imbalance is exactly same as the classical 
equation of motion of a particle at q in GL po- 
tential Vgl (q) [13 . This GL potential changes its 
shape when E/S > — , and two minima appear 
for non vanishing population imbalance. For MST, 
the classical particle has total energy less than the 
barrier height in order to be localized at one of the 
minima of the GL potential. Two minima of the GL 



potential at population imbalance q = ± y 1 — 

correspond to the fixed points of MST states with 
energy E/S = h(a + 1/a). The frequency of har- 
monic oscillation around the potential minimum is 
given by lj/J = y/cr — 1. For a > 1, the GL po- 
tential can have two minima around which dynam- 
ically stable MST state can exist. 

From the above classical analysis we notice that a c = 1 
is a critical coupling strength which separates two types 
of dynamical states. Next we analyze the dynamical 
and spectral properties of the effective spin Hamiltonian 
(Eq. @ ) quantum mechanically. In order to analyze var- 
ious dynamical branches discussed above, we study nu- 
merically full quantum dynamics of BJJ with finite (but 
large) number of atoms and compare it with classical dy- 
namics. We use the effective spin Hamiltonian (in Eq.®) 
to describe the quantum dynamics of BJJ. Any quantum 
mechanical state of the system can be written as a linear 
combination of 2S + 1 basis states, 
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where \a) is the eigensate of the operator S z with eigen- 
value a, and a a s are complex (in general time depen- 
dent) amplitude. To study the quantum dynamics, we 



FIG. 1: Oscillation of fractional population imbalance AN/N 
with time t (in units of 1/J) for different values of coupling 
constant a and initial conditions. Results obtained from 
quantum dynamics(solid lines) are compared with classical 
dynamics (dotted line) for same initial conditions, a) Ac 
Josephson oscillation for N = 100, a = 0.5, and with clas- 
sical initial condition cos 9 = 0.2, <j> = 0. b)7r-Oscillation for 
N = 100, a = 0.5, cos# = 0.2, and = vr. c) Self trapping 
for TV = 100. a = 1.5. cos (9 = 0.85. and d> = n. d) DamDins 
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time evolve any initial state by the time-evolution op- 
erator e~ ttH , and then evaluate the expectation values 
of the observablcs. To compare the quantum dynamics 
with the classical dynamics of the spin, we take the ini- 
tial state of the spin as a coherent state corresponding to 
its classical orientation . Classical nature of a large spin 
is well described by a spin coherent statefl8j. 
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1 



E 



(10) 



where, the complex number £ = tan represents the 
orientation of the corresponding classical spin vector. We 
expand the initial state in terms of the eigen-vectors of 
the Hamiltonian and obtain the final state at time t, mul- 
tiplying each component by a factor e %Ent where E n is 
corresponding eigenvalue. For N = 100 (S = N/2), we 
calculate the time evolution of fractional population im- 
balance (Nl — N2)/N between two condensates (with 
average number of particles Nl and N2) in BJJ, from 
the expectation value of the operator S z /S, and com- 
pare it with the classical result obtained from time inte- 
gration of Eqn.dU and Eqn.JS]). To study different types 
of dynamics, we choose the initial values of the classical 
variables close to the fixed points and the initial coher- 
ent state can be obtained from the corresponding value 
of £. Three different types of Joscphson dynamics (as 
discussed above) are shown in Fig. 2, for different values 
of the coupling constant a and are compared with the 
classical dynamics. The ac Joscphson dynamics and V- 
oscillation' for N = 100, show very good agreement with 
the classical dynamics. However the quantum dynamics 
around the MST state shows some deviation from the 
classical dynamics, also it shows some beating effect (see 
fig.lc). For the same initial number-imbalance but with 
phase (j) = 0, quantum dynamics shows a damping and 
long time revival phenomena as shown in Figld. 



III. EFFECTIVE POTENTIAL METHOD AND 
MACROSCOPIC QUANTUM TUNNELING 

In this section we consider interesting aspects related 
to the collective quantum effects in a BJJ in terms of the 
'effective potential method' of a large spin[l9j]. The 'ef- 
fective potential' of the quantum spin Hamiltonian gives 
a clear picture about various (meta) stable dynamical 
branches and also elucidate possible collective quantum 
tunneling in a BJJ. 

First we outline the method of effective potential for 
the spin Hamiltonian in Eq.@ describing the BJJ. The 
eigenvalue equation of the Hamiltonian in Eq.([2]) can be 
written as, 



Ea 



U 



<y 2 a n 



J 



2 [^(S + a)(S- a +1K_! 
<j+l)a a+1 ], (11) 




FIG. 2: Effective potential: V(x) in the range — ir < x < n, 
for a) a = 0.4, and b) a = 1.4. W(x) for c)q = 0.4 and d) 
a = 1.4. 



where the eigenstates are expanded in terms of the basis 
states of S z , as described in Eq.©. By introducing a 



generating function $(x) 



of an 



^/(S+crV.(S-a)\ 

auxiliary variable x, above eigenvalue equation Eq. (|lip 
can be written as a differential equation [l9(, 



Ud 2 <S> J 



2 dx 2 2 



2S cosx$(x) — 2sinx 



<9$ 

dx 



(12) 



It is interesting to note that $(x) = (£|V>)e lS ^, with 
£ = e ix . So the variable x can be interpreted as the phase 
4> between two condensates with equal number of parti- 
cles in each trap i.e |£| = 1 (or cos 9 = 0). The transfor- 



mation $(.t) = e 



-J cos x/U 



x{ x )i reduces the above eigen- 



value equation to an effective Schrbdingcr's equation, 



a 



2(s + l/2) 



;d 2 xX + V(x) X 



(13) 



where, £ = E/(J(S+l/2)), a' = U(S + l/2)/J « a, and 
the effective potential is given by, 



1 



V(x) = sin x — cos a;. 

2a' 



(14) 



First 25*+ 1 eigenvalues of the above Schrbdinger's equa- 
tion are the eigenvalues of the original spin Hamiltonian 
Eq.([2]). Above effective potential V(x) for the relative 
phase can also be derived from the Eq. ((H) and Eq.© de- 
scribing the classical dynamics. The relative phase <f> 
between the condensates obeys the Newton's equation of 
motion in the effective potential V(<fi) [l2[. Above Hamil- 
tonian in Eq. (fl~3"| with the effective potential Eq. (fT4")) is 
an exact quantum phase description (or 'rotor model') of 
Josephson junctions [20|. 

It is interesting to note following significant points of 
the above equation: 

• The effective Planck constant in Eq. (fT"3]) is 1/5. 
Which implies that quantum fluctuation is sup- 
pressed for large spin systems (or large number of 
particles in the trap). 
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• The role of effective mass is played by the inverse 
of the coupling constant 1/a. With increasing cou- 
pling constant a the fluctuation in number imbal- 
ance (or S z ) is reduced and the quantum fluctua- 
tion in phase </> is enhanced. 

• As depicted in Fig. 2a and Fig. 2b, the shape of the 
effective potential V(x) changes at a critical cou- 
pling a = 1. For a < 1, the potential has two 
minima at x = 2nir and x = (2n + 1)tt. These 
two minima correspond to the stable phase oscilla- 
tions between the condensates discussed in previous 
section. The 7r-oscillation corresponds to the local 
minima of the effective potential at x = (2n + l)ir. 
The phase fluctuations around the classical steady 
state solutions can be obtained from harmonic ap- 
proximation of the potential around the local min- 
ima. The potential at local minima x± can be ap- 
proximated as, 



1 



2a 



(15) 



where, x± corresponds to the points and ±7r re- 
spectively. Note that the frequency of the harmonic 
potential exactly matches with the frequency of 
small amplitude oscillations. The phase fluctua- 
tions is given by, 



({x-x ± f) 



2SVT±~a~ 



(16) 



Near the critical point the phase-fluctuation of ir- 
oscillation becomes very large. 



A. Decay of the V-oscillation' 

As it is evident from the effective potential V(x), the 
state with a relative phase difference it between the con- 
densates is a mctastable state which is separated from 
the ground state by a potential barrier (see Fig. 2a). A 
metastable state of many particles can decay by means of 
collective quantum tunneling through the barrier, known 
as 'macroscopic quantum tunneling' fMQT) (2~i| . Since 
this collective state of many bosons is represented by a 
point particle in a mcta-stable local minima of the effec- 
tive potential V(x), we can use simple quantum mechan- 
ical tunneling formula to calculate the collective decay 
rate. The tunneling probability of a single quantum par- 
ticle through a potential barrier scales as ~ e ~\f-^r*- L / h ^ 
where V max is the height of the potential barrier, L is the 
length of the barrier and m is the mass of the particle. 
In case of quantum phase model described by Eq. (flTJ|) . 
the effective Planck constant h scales by 1/5 and hence 
the collective tunneling rate is suppressed by a factor 
of S (or by the total number of particles A)in the ex- 
ponential. As seen from the simple scaling analysis, the 
decay of the metastable '7r-oscillation' state is suppressed 



exponentially by the number of particles N, and is neg- 
ligible except in the region close to the critical coupling 
strength a ~ 1 where the barrier height becomes vanish- 
ingly small. The potential near the meta stable state of 
±7r-phase can be approximated as, 



e C.2 1 ~4 



V{x) = 1 + -x 2 --x 



(17) 



where x = x — x—, and e = 1 — a can be considered as 
a small parameter close to the critical point. To calcu- 
late the imaginary part of the ground state energy of the 
metastable state, we use the imaginary time path integral 
method. The imaginary time classical path in the quartic 
potcntial(in Eq. (|17[) ) is known as instanton solution, 



x(t) = 2i/e sech [y/e(r - To)] 



(18) 



where r denotes imaginary time. Above equation de- 
scribes the classical path of a particle in the inverted 
potential given in Eq. ([T7)) . where the particle starts from 
x = X— at r = — oo, reaches the turning point x = 2y/e at 
t = tq, and returns to the original point at r = oo. The 
imaginary time 'action' S c corresponding to this classical 
path is given by, 



S c = 



16 



(It 



3/2 



(^r+v(x)-v(x.) 



(19) 



By using instanton technique, the decay rate T of a quan- 
tum particle in an unstable quartic potential ( given 
by Eq. (fTT)) ) can be calculated in terms of the imagi- 
nary action S c and the frequency of oscillation u> at 
the metastable local minimum, and the analytical ex- 
pression is T = ujhy/6S^7rhe~ Sa ^ h ^^ l . In the 'effec- 
tive' Schrddinger equation approach of the spin model 
as written in Eq . (fT3|) , h is replaced by 1/S, and the 
frequency of oscillation at mctastable state x = ±tt is 
uj = yl — a m yfe. Substituting these parameters in 
the above expression, we obtain the decay rate of the 
metastable 7r-oscillation close to the critical point (when 
e« 1), 



96 
3irS 



; 5/2 p -16Se 3 / 2 /3 



(20) 



B. Energy spectrum 

Similar to the change in shape of the potential V(x), 
we also notice a distinct change in the energy spectrum 
of the spin Hamiltonian in Eq.([2]), as the coupling con- 
stant a is varied across the critical value a c = 1. We also 
analyze the energy spectrum of the Hamiltonian given 
in Eq.([5]) by numerical diagonalization, which is shown 
in the Fig. 3 ( a and b) for S = 50. The interesting 
feature in the energy spectrum is appearance of doubly 
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FIG. 3: a)Energy eigenvalues of the Hamiltonian in Eq.© 
for N = 50, a = 0.4(circles connected with line) and a — 
2.5(squares). b)Energy gap AE/S as a function of coupling 
strength a for N = 500 (solid line) and N = 800 (dotted line). 



(quasi) degenerate excited states above the critical cou- 
pling strength a > 1 (see Fig. 3(b)). First the degeneracy 
appears in the highest energy level, and then with the in- 
creasing coupling strength a, other excited states become 
pairwise degenerate. To quantify this transition in the 
energy spectrum, we calculate the energy gap A between 
the highest energy level E2S+1, and the next one E2S 
as a function of the coupling constant a (as depicted in 
Fig. 2(c)). For finite S, the energy gap A vanishes slightly 
above the critical coupling a = 1. With the increasing 
values of S, the transition becomes sharper and the gap 
vanishes closer to the critical coupling. This transition 
can be explained in terms of the 'effective potential' ap- 
proach of the spin system and is directly linked to the 
'self-trapping' phenomena. 

The energy spectrum of the original spin Hamiltonian 
in Eq.([2]) is exactly same as that of the Schrodinger's 
equation with an effective potential given in Eq.(|13j). and 
Eg. (1141) . By performing two successive transformations, 
x — x' + 7r, and x' = ty, we obtain following eigenvalue 
equation from Eq. (fT3|) . 



2(s + l/2) 



2 d 2 v x(v) + w(y)x(y), (21) 



where, the effective potential W(y) = j- sinh 2 y — coshy, 
and k = —£. Above transformation leads to the mapping 
of the highest eigenvalue -E2S+1 of the Sclrrddinger equa- 
tion given in Eq. (|T3"|) to the ground state of the effective 
potential W(y) in Eq. (j2"T1) . For a < 1, the effective po- 
tential has only one minimum and then it takes a shape 
of a double well for a > 1, as shown in Fig. 2(c) and 
Fig. 2(d). Double well structure of the effective potential 
W explains the appearance of doubly degenerate excited 
states of the original spin Hamiltonian for the coupling 
strength a > 1. When the kinetic energy term in the 
Hamiltonian (Eq. (j2"Tj) is neglected, the classical ground 
state energy can be obtained from the potential W at 
two degenerate minima. Classical ground state energy 
h(a + — ) of the double well potential is exactly same as 
that of the fixed point corresponding to the MST state 
with population imbalance ±s/T— 1/a 2 . As discussed 
earlier, the MST state can be described by the dynam- 
ics of a fictitious classical particle in the GL potential 
(in Eq.((5])) with energy E less than the barrier height. 



The formation of the doubly degenerate excited states 
in the quantum mechanical spectrum can be understood 
from the periodic orbits of the classical motion near two 
minima of the GL potential. 



Energy splitting and Macroscopic quantum 
tunneling between MST states 



For S — > 00, the kinetic energy term in Eq. (|21j) van- 
ishes and two energy sates corresponding to the min- 
ima of the potential W(y) are degenerate. But for fi- 
nite values of S, the fictitious quantum particle described 
by Schrodinger's equation in Eq. ([21[) , can tunnel be- 
tween the minima of the potential and produce a small 
splitting between the highest degenerate eigenvalues dis- 
cussed above. The tunneling time T of the particle can 
be estimated from the relation T ~ h/AE, where AE 
is the energy splitting. This simple picture elucidates 
a novel phenomenon of MQT of the condensate in BJJ 
between the MST with equal classical energy. 

Next, we study the MQT of the condensate from one 

fixed point of MST with number imbalance */l — 

to the other fixed point \f \ — I /a 2 , which are repre- 
sented by the minima of 'effective potential' W. The 
ground state energy splitting due to the quantum tun- 
neling can be calculated by the instanton technique. 
The tunnel splitting of the ground state energy of the 
Schrodinger's equation in Eq.([2"Tj) has already been cal- 
culated by several authors in the context of uni-axial 
paramagnet [H, [23|. For a > 1, the effective po- 
tential develops two minima which satisfy the equation 
cosh(a; TO j n ) = a. The imaginary time classical path join- 
ing two minima of the potential can be calculated exactly 



tanh- 



1 , 1 
— tanh— lot, 
1 2 ' 



(22) 



where to = \/a 2 — 1 is the frequency of oscillation at the 
minima of the potential. Classical action corresponding 
to this instanton solution is given by, 



S c = 2ln a + vV - 1 



2VV 



1/a (23) 



From the well known semiclassical formula for tunnel 
splitting[24],[25j], we obtain the energy difference between 
two MST states corresponding to the minima of the ef- 
fective potential W, 



AE 



1) 



[a + Va 2 - lj 



2S 



D 2SVa 2 -l/a 



(24) 



Apart from the non-linear oscillations around the self 
trapped state, the BJJ can show a very interesting quan- 
tum oscillation from one MST state to other MST state 
with equal and opposite number imbalance, and the es- 
timated time period of such oscillation due to MQT is 



7 



~ sSb- For large number of particle (or large S), the 
tunneling time will be very large and the system will 
practically be trapped to one of the MST states depend- 
ing on the initial condition. But for finite number of par- 
ticle in BJJ, and by tuning the coupling constant close 
to the critical point, such MQT can be observed. 

IV. SUMMARY 

To summarize, we consider two weakly coupled con- 
densates described by an effective two-site BHM and 
studied quantum dynamics and collective quantum ef- 
fects. We study quantum mechanically various dynami- 
cal branches of the Joscphson oscillations and transition 
between them. Two-site BHM with N bosons can be 
mapped on to a simple model describing a quantum par- 
ticle in an 'effective potential' where the Planck constant 
h is reduced by a factor of N/2. This 'effective poten- 
tial' gives a clear picture of the '7r-oscillation' and macro- 
scopic self-trapping phenomena in BJJ. At the transition 
between these two dynamical branches, the 'effective po- 
tential' changes its shape. 

Apart from the dynamics, the novel quantum effects 
like collective decay and MQT emerge from the simple 
picture of 'effective potential' method. The '7r-oscillation' 
of the condensate corresponds to a meta stable minimum 
of the potential V(x), and can decay by means of MQT 



through the potential barrier. Scmiclassically we calcu- 
late the decay rate T for small potential barrier and ob- 
tain a power law behavior T ~ (a c ~ a) 5 / 4 close to the 
critical coupling strength a c . 

We also notice a remarkable feature in the energy spec- 
trum as the coupling constant a changes across the crit- 
ical value. Quasi-degenerate pairs of excited states ap- 
pear in the energy spectrum for coupling strength larger 
than the critical value a c . The doubly degenerate high- 
est energy states can be interpreted in terms of an effec- 
tive double well potential W{x) (shown in Fig. 2(d)), and 
they correspond to the classically degenerate MST states 
with equal and opposite number imbalance. We study 
the MQT between these two MST states and estimate 
the tunneling time from the energy splitting. The MQT 
rate is exponentially suppressed by number of particles 
N, nevertheless these novel macroscopic quantum phe- 
nomena can be observed in a BJJ with finite number of 
atoms and by tuning the coupling constant close to the 
critical value. 
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